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Abstract 

We illustrate the effects of one-loop weak corrections onto the production of neutral 
gauge bosons of the Standard Model at RHIC-Spin, Tevatron and LHC, in presence of 
quark/gluon radiation from the initial state. We find such effects to be rather large, up 
to O(10 — 20%) in typical observables at all such colliders, where the cross section is 
measurable, thus advocating their inclusion in precision analyses. 
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1 Prompt-photon and on-shell Z hadro-production 



The neutral-current processes (V = 7, Z) 

qq -> gV and q(q)g -> g(g)V (1) 

with V" — > are two of the cleanest probes of the partonic content of (anti)protons, in 

particular of antiquark and gluon densities. In order to measure the latter it is necessary 
to study the vector boson p T spectrum. According to [1, 2] the gluon density dominates for 
Pt > Q/2 where Q is the lepton pair invariant mass. In the presence of polarised beams 
these reactions give access to the spin-dependent gluon distribution which is presently only 
poorly known. Thanks to the introduction of improved algorithms [3]- [5] for the selection of 
(prompt) photons generated in the hard scatterings (1), as opposed to those generated in the 
fragmentation of the accompanying gluon/quark jet, and to the high experimental resolution 
achievable in reconstructing Z — > (£ = e, //) decays, they are regarded - together with 
the twin charged- current channels 

qq'^gW and q(q)g - q'ffiW, (2) 

wherein W — > Iv? - as precision observables in hadronic physics. In fact, in some instances, 
accuracies of order one percent are expected to be attained in measuring these processes [6], 
both at present and future proton-(anti)proton experiments. These include the Relativistic 
Heavy Ion Collider running with polarised proton beams (RHIC-Spin) at BNL {\fs pp = 300 — 
600 GeV), the Tevatron collider at FNAL (Run 2, y/s pp = 2 TeV) and the Large Hadron 
Collider (LHC) at CERN (^7s pp = 14 TeV). 

Not surprisingly then, a lot of effort has been spent over the years in computing higher 
order corrections to all such Drell-Yan type processes. To stay with the neutral-current ones, 
the subject of this paper, these include next-to-leading order (NLO) QCD calculations of both 
prompt-photon [7, 8] and vector boson production [9]. QCD corrections to the pt distributions 
have been computed in Refs. [10, 11]. As for the full 0(a) Electro- Weak (EW) corrections to 
Z production and continuum neutral-current processes (at zero transverse momentum), these 
have been completed in [12] (see also [13]), building on the calculation of the QED part in 
[14]. 

In the case of polarised (anti)proton beams, the process of calculating higher order correc- 
tions has proceeded more slowly [15, 16]. NLO QCD corrections to the transverse momentum 
spectrum of Drell-Yan type processes via neutral-currents in presence of (longitudinal) spin ef- 
fects from the initial state can be found for the non-singlet case in [17, 18], while the complete 
calculation has been recently published in Ref. [19] (see also Ref. [20]). 

The relatively large impact of one-loop EW corrections, as compared to the QCD ones, can 
be understood (see Refs. [21]-[22] and references therein for reviews) in terms of the so-called 
Sudakov (leading) logarithms of the form aw log 2 (v / 5/M^), which appear in the presence of 
higher order weak corrections (hereafter, a w = oem/ sin26* w , with a EM the Electro-Magnetic 
(EM) coupling constant and # w the weak mixing angle) 2 . These 'double logs' are due to a 

In some cases, leading (~ ofy log 2 " (s/M w 2)), sub-leading (~ d^\og 2n ~ l {s/M w 2)) and sub-sub-leading 
(~ oeyj \og 2n ~ 2 (s / M w 2)) can be resummed. 
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lack of cancellation of infrared (both soft and collinear) virtual and real emission in higher 
order contributions due to lU-exchange in spontaneously broken non-Abelian theories. 

The problem is, in principle, present also in QCD. In practice, however, it has no observ- 
able consequences, because of the averaging on the colour degrees of freedom of partons, forced 
by their confinement into colourless hadrons. This does not occur in the EW case, where, 
e.g., the initial state can have a non-Abelian charge, dictated by the given collider beam 
configuration. Modulo the effects of the Parton Distribution Functions (PDFs), which spoil 
the subtle cancellations among subprocesses with opposite non-Abelian charge, for example, 
this argument holds for an initial quark doublet in proton-(anti)proton scatterings. These 
logarithmic corrections (unless the EW process is mass-suppressed) are universal (i.e., process 
independent) and are finite (unlike in QCD), as the masses of the EW gauge bosons provide 
a physical cut-off for lU-boson emission. Hence, for typical experimental resolutions, softly 
and collinearly emitted weak bosons need not be included in the production cross-section and 
one can restrict oneself to the calculation of weak effects originating from virtual corrections. 
In fact, one should recall that real weak bosons are unstable and decay into high transverse 
momentum leptons and/or jets, which are normally captured by the detectors. In the def- 
inition of an exclusive cross section then, one tends to remove events with such additional 
particles. Under such circumstances, the (virtual) exchange of Z-bosons also generates simi- 
lar logarithmic corrections, a^flog (s/l/M^). Besides, the genuinely weak contributions can 
be isolated in a gauge-invariant manner from purely EM effects, at least in some simpler cases 
- which do include processes (1) but not (2) - and the latter may or may not be included in 
the calculation, depending on the observable being studied. 

A further aspect that should be recalled is that weak corrections naturally introduce parity- 
violating effects in observables, detectable through asymmetries in the cross-section, which 
are often regarded as an indication of physics beyond the Standard Model (SM) [6, 23, 24]. 
These effects are further enhanced if polarisation of the incoming beams is exploited, such 
as at RHIC-Spin [25, 26]. Comparison of theoretical predictions involving parity- violation 
with experimental data is thus used as another powerful tool for confirming or disproving 
the existence of some beyond the SM scenarios, such as those involving right-handed weak 
currents [27], contact interactions [28] and/or new massive gauge bosons [29, 30]. 

In view of all this, it becomes of crucial importance to assess the quantitative relevance 
of weak corrections affecting processes (l)-(2). It is the aim of our paper to report on the 
computation of the full one-loop weak effects entering processes (1) while the study of those 
for (2) will be deferred to a future publication [31]. 

2 Calculation and results 

Since we are considering weak corrections that may be identified via their induced parity- 
violating effects and since we wish to apply our results to the case of polarised proton beams, it 
is convenient to work in terms of helicity Matrix Elements (MEs). Here, we define the helicity 
amplitudes by using the formalism discussed in Ref. [32]. At one-loop level such helicity 
amplitudes acquire higher order corrections from: (i) self-energy insertions on the fermions 
and gauge bosons; (ii) vertex corrections and (iii) box diagrams. All such contributions are 
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pictured in Fig. 1. 

The self-energy and vertex correction graphs contain ultraviolet divergences that have 
been subtracted here by using the 'modified' Minimal Subtraction (MS) scheme at the scale 
\i = M Z - Thus the couplings are taken to be those relevant for such a subtraction: e.g., 
the EM coupling, qem, has been taken to be 1/128 at the above subtraction point. The one 
exception to this renormalisation scheme has been the case of the self-energy insertions on 
external fermion lines, which have been subtracted on mass-shell, so that the external fermion 
fields create or destroy particle states with the correct normalisation. 

All these graphs are infrared and collinear convergent so that they may be expressed 
in terms of Passarino-Veltman [33] functions which are then evaluated numerically. The 
expressions for each of these diagrams have been calculated using FORM [34] and checked 
by an independent program based on FeynCalc [35]. For the numerical evaluation of the 
scalar integrals we have relied on the FORTRAN package FF [36]. A further check on our 
results has been carried out by setting the polarisation vector of the V^-boson proportional to 
its momentum and verifying that the sum of all one-loop diagrams vanishes, as required by 
gauge and BRST invariance. The full expressions for the contributions from these graphs are 
too lengthy to be reproduced here. 

In both processes in (1), external (anti)quarks have been taken massless and both vector 
bosons (V = 7, Z) have been put on-shell. In contrast, the top quark entering the loops in 
both reactions has been assumed to have the mass m t = 175 GeV. The Z mass used was 
M z = 91.19 GeV and was related to the I^-mass, M w , via the SM formula M w = M z cos9 w , 
where sin 29 w = 0.232. (Corresponding widths were T z = 2.5 GeV and T w = 2.08 GeV.) 
For the strong coupling constant, as, we have used the one-loop expression at Q 2 = \fl with 
A^- -4 '* chosen to match the value required by the (LO) PDFs used. The latter were Gehrmann- 
Stirling set A for RHIC and Martin-Roberts-Stirling-Thorne set 2001 LO for Tevatron/LHC 
[37]. 

2.1 RHIC 

The following beam asymmetries can, e.g., be defined at RHIC-Spin: 

All da = da ++ — da^ 

+ da — da h , 

A L da = da_ — da + , 

Apy da = da — da ++ . (3) 

The first two are parity-conserving while the last two are parity-violating. Figs. 2-3 show the 
NLO distributions in such quantities for both 7 and Z final states, alongside those for the total 
cross section, as a function of the transverse momentum p?, within the pseudo-rapidity range 
\y\ < 1. The corrections due to full one-loop weak effects are also presented (for the cases 
in which the Born level result is non-zero). Effects onto the total cross sections are rather 
small, below the percent level, as expected, because of the rather low centre-of-mass (CM) 
energy available at partonic level, which is comparable with M w and M z , so that logarithmic 
corrections are not enhanced. Nonetheless, in the case of 7-production, at p T = 10 GeV, the 
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total NLO yearly rates of approximately 150,000 and 1,350,000 events accessible at low and 
high energy, respectively, contain a sizable contribution due to purely weak effects, about 200 
and 1,500 events in correspondence of y/s = 300 and 600 GeV (for the values of luminosity 
200 and 800 pb _1 , respectively). In case of Z-pro duct ion, only at 600 GeV NLO effects are 
sizable, as they are responsible for 3 events being subtracted (the correction is negative) to 
the LO prediction of 457 events (at pt = 10 GeV). (The LO rate at 300 GeV for Z-production 
at such a transverse momentum is of only 7 events, unaffected by NLO weak effects.) 

Rather large effects do appear in general for the asymmetries, particularly for the case of 
Z boson final states. In fact, for the latter, in the case of Al and All, they can range up to 
±(15 — 20)%, while they are somewhat lower for the case of Apy, 5% or so. Unfortunately, 
none of such NLO effects on the asymmetries is detectable, because of the poor production 
rate of Z-bosons. For photonic final states, one-loop weak effects on A LL are not much larger 
than those on the total rates, nonetheless, they might just be observable at low p T . For the 
cases of A L and A PV , which for the photon are exactly zero at Born level in massless QCD, 
one-loop weak effects are too poor to be observed experimentally. 

2.2 Tevatron and LHC 

Figs. 6-5 show the effects of the 0(asa^ w ) terms relatively to the C(q;sq;ew) Born results 
(«em replaces «ew for photons), as well as the absolute magnitude of the latter, as a function 
of the transverse momentum, at Tevatron and LHC, respectively. The corrections are found 
to be rather large at both colliders, particularly for Z-pro duct ion. In case of the latter, such 
effects are of order -7% at Tevatron for p T « 300 GeV and -14% at LHC for p T « 500 GeV. 
In general, above px ~ 100 GeV, they tend to (negatively) increase, more or less linearly, 
with pt- Such effects will be hard to observe at Tevatron but will indeed be observable at 
LHC. For example, at FNAL, for Z-production and decay into electrons and muons with 
BR(Z —> e, fi) 6.5%, assuming L = 2 — 20 fb _1 as integrated luminosity, in a window of 10 
GeV at p T = 100 GeV, one finds 500-5000 Z + j events at LO, hence a 5a /a « -1.2% EW 
NLO correction corresponds to only 6-60 fewer events. At CERN, for the same production 
and decay channel, assuming now L = 30 fb" 1 , in a window of 40 GeV at pt = 450 GeV, 
we expect about 2000 Z + j events from LO, so that a 5a j a —12% EW NLO correction 
corresponds to 240 fewer events. In line with the normalisations seen in the top frames of 
Figs. 6-5 and the size of the corrections in the bottom ones, absolute rates for the photon are 
similar to those for the massive gauge boson while 0(asa^ w ) corrections are about a factor 
of two smaller. 

3 Conclusions 

Altogether, the results presented here point to the relevance of one-loop O^a^a^j) weak contri- 
butions for precision analyses of prompt-photon and neutral Drell-Yan events at both Tevatron 
and LHC, also recalling that the residual scale dependence of the known higher order QCD 
corrections to processes of the type (1) is very small in comparison [11]. Another relevant 
aspect is that such higher order weak terms introduce parity-violating effects in hadronic ob- 
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servables [26], which might just be observable at (polarised) RHIC-Spin, particularly in the 
case of photons. (The case for polarised beams at the LHC, enabling the study of parity- 
violating asymmetries on the same footing as at RHIC-Spin, is currently being discussed as 
one of the possible upgrades of the CERN collider.) The size of the mentioned corrections, 
relative to the lowest order results, is insensitive to the choice of PDFs. EM effects were 
neglected here because they are not subject to logarithmic enhancement, thus smaller with 
respect to the weak ones, or parity-violating effects either. However, their computation is 
currently in progress [31]. 
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Figure 1: Graphs describing processes (1) in the presence of one-loop weak corrections. The shaded 
blob represents all the contributions to the gauge boson self-energy and is dependent on the Higgs 
mass (we have set Mh = 115 GeV). We neglect loops involving the Higgs boson coupling to the 
fermion line. The graphs in which the exchanged gauge boson is a VF-boson are accompanied by 
those in which the latter is replaced by its corresponding Goldstone boson. There is a similar set of 
diagrams in which the direction of the fermion line is reversed, with the exception of the last graph, 
as here reversal does not lead to a distinct topology. 
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Figure 2: The transverse momentum dependence of the 7-boson cross section in (1) as well as of 
the beam asymmetries in (3) at NLO (large frames) and the size of the one-loop weak corrections 
(small frames, limitedly to case in which the latter is non-zero), at RHIC-Spin {yfs pp = 300 and 600 
GeV). Notice that the pseudorapidity range of either particle in the final state is limited to \rj\ < 1. 



9 



> 

0) 

O 

\ 

E-i 
ft 



b 



solid 



10 L 



Z-production (rj < 1) 
6 = (NL0-L0)/L0 
V s = 300 GeV (NLO) 



dashed 

i i i 



V s = 600 GeV (NLO) 



10 
10" 
10' 



-3 
-6 
-9 




0.0 
-0.2 
-0.4 
-0.6 
-0.8 

i 

0.125 
0.100 
0.075 
0.050 
0.025 
0.000 

10 

5 

-5 
-10 
-15 



= 1111 


1 1 1 1 


II 


II 


£ 


= 1111 


1 1 1 1 


II 


Mill 1 nN 


. 1 t 



20 
10 


-10 
-20 



-III 


II 


1 1 1 ^ II II 


III 


Illllllt 


=T1 1 


II 


If 

V--7 
1 1 1 II II 


III 





20 40 



60 80 100 

Pt (GeV) 







20 40 



-Mil 


1 1 1 1 


ii 


ii 


mi 






-i-rfi 


ii 


ii 


Ml 





= 1111 


MJM 1 


1 1 1 | 1 1 1 


III 


t 




/ 








= \ 
= 1 1 M 


1 
1 

1 

j I M I 


M I 1 M I 


Ml 


f 







20 40 



60 80 100 
Pt (GeV) 



0.125 
0.100 
0.075 
0.050 
0.025 
0.000 

10.0 

7.5 
5.0 
2.5 
0.0 
-2.5 
-5.0 



60 80 
Pt (GeV) 



100 



-Ml 


1 II 1 1 1 


mi 


III 


A- 




ri+rn ii 


Ml 


III 





> 

13 



20 40 



60 80 
Pt (GeV) 



100 



Figure 3: The transverse momentum dependence of the Z-boson cross section in (1) as well as of 
the beam asymmetries in (3) at NLO (large frames) and the size of the one-loop weak corrections 
(small frames), at RHIC-Spin {\fs pp = 300 and 600 GeV). Notice that the pseudorapidity range of 
either particle in the final state is limited to \rj\ < 1. 
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Figure 4: The transverse momentum dependence of the 7- and Z-boson cross sections in (1) at LO 
(top frame) and the size of the one-loop weak corrections (bottom frame), at Tevatron (y/s p = 2 
TeV). Notice that the pseudorapidity range of the jet in the final state is limited to \rj\ < 3. 



11 



I Pt LHCLO I 




I 1 

200 400 600 800 1000 1200 1400 1600 1800 
Pt (GeV) 



| LHC (NLO-LO)/LO | 
0.02 1 




U I I I I I I ! I I ! I I I I ! I I I 

50 100 150 200 250 300 350 400 450 500 

Pt (GeV) 



Figure 5: The transverse momentum dependence of the 7- and Z-boson cross sections in (1) at LO 
(top frame) and the size of the one-loop weak corrections (bottom frame), at LHC {\fs pp = 14 TeV). 
Notice that the pseudorapidity range of the jet in the final state is limited to \rj\ < 4.5. 
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Erratum to: 
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The results presented in Fig. 4 of the original paper mistakenly refer to a pp collider of 
\fs pp = 2 TeV instead of a pp one. The correct results for the effects of the O(aso^ w ) 
terms relatively to the O^asa^w) Born results («em replaces «ew f° r photons), as well as the 
absolute magnitude of the latter, as a function of the transverse momentum at Tevatron are 
shown in Fig. 6 below. The corrections are of order -6% for Z + j production at Tevatron 
for px ~ 300 GeV. Since both the size of the corrections and the cross section for moderate 
values of pt are similar to those for a pp collider, our conclusions that such effects will be hard 
to observe at Tevatron but will indeed be observable at LHC are unchanged. 
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Figure 6: The transverse momentum dependence of the 7- and Z-boson cross sections in qq — * 
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